Nonlinear Plasma Wave in Magnetized Plasmas 
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Abstract Nonlinear axisymmetric cylindri- 
cal plasma oscillations in magnetized collision- 
less plasmas are a model for the electron fluid 
collapse on the axis behind an ultrashort rela- 
tivisically intense laser pulse exciting a plasma 
wake wave. We present an analytical descrip- 
tion of the strongly nonlinear oscillations show- 
ing that the magnetic field prevents closing of 
the cavity formed behind the laser pulse. This 
effect is demonstrated with 3D PIC simulations 
of the laser-plasma interaction. An analysis of 
the betatron oscillations of fast electrons in the 
presence of the magnetic field reveals a charac- 
teristic "Four- Ray Star" pattern which has been 
observed in the image of the electron bunch 
in experiments [T. Hosokai, et a/., Phys. Rev. 
Lett. 97, 075004 (2006)]. 



I. INTRODUCTION 

Theoretical and experimental studies of elec- 
tromagnetic wave propagation in magnetized 
plasmas are of key importance for a vast range 
of problems in space and laboratory physics 
[U 12]. The interest towards relativistically 
strong plasma waves stems both from astro- 
physical applications and from high field sci- 
ence, i. e. from relativistic plasmas research 
[3]. Relativistic plasma waves can be gener- 
ated as wake waves in a plasma behind an ul- 
trashort laser pulse [11 [5] or by bunches of ultra- 
relativistic electrons [6^, providing extremely 



strong electric fields used for the acceleration of 
high energy charged particles. Electron acceler- 
ation by the plasma wake field above energies of 
GeV-ies has been demonstrated in a number of 
experiments [7 . Physical processes associated 
with the nonlinear wake wave play a key role in 
high-harmonic generation [8 and in many other 
aspects of laser-plasma physics [3 , including the 
frequency up-shifting and intensification of the 
electromagnetic radiation [9l [10] . 
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An important property of a nonlinear wake 
wave is that when its amplitude exceeds the 
wavebreaking threshold the wake breaks trans- 
ferring its energy to the plasma electrons, thus, 
on one hand, providing a robust mechanism of 
electron injection to the acceleration phase of 
the wake field [11 and forming regular nonlinear 
structures moving with relativistic velocity [10] 
and, on the other hand, imposing a constraint 
on the achievable wake wave amplitude [12(113]. 
A tightly focused and sufficiently intense laser 
pulse excites a wake wave in the regime, where 
in addition to a longitudinal push the laser pulse 
expels electrons also in transverse direction [13]- 
16 , forming a cavity void of electrons in the 
first period of the wake wave. The parameters 
of the cavity in the electron density behind the 
laser pulse can be deduced from the fact that 
all the electrons for a sufficiently intense laser 
pulse are pushed aside resulting in a multiffow 
motion owing to the collisionless nature of the 
low-density plasma [TSJ [16] . The region void of 
electrons is positively charged and attracts un- 
perturbed electrons on the side. In the near axis 
region, where the cavity is closing, the electron 
fluid dynamics can be modeled by that of an 
axisymmetric cylindrical Langmuir wave. Non- 



linear cylindrical and spherical Langmuir waves 
have been considered in Refs. [17], where it 
was shown that the finite amplitude electro- 
static waves break due to nonlinear effects and 
the non-planar geometry. 

Below we report the results of analytical in- 
vestigations and computer simulations of non- 
linear Langmuir oscillations in a plasma with 
a relatively strong homogeneous axial magnetic 
field. There are several mechanisms of strong 
magnetic field generation in the laser matter in- 
teraction. Large amplitude magnetic fields are 
generated by temperature gradients near the 
laser focal spot, by the hot electron currents 
driven by an ultraintense laser pulse, in the 
plasma jets induced by them, and by the com- 
pression of thin shells by the laser pulses fTS^. 
Axial magnetic fields can be generated by cir- 
cularly polarized laser radiation via the Inverse 
Faraday Effect [19^. Strong magnetic fields can 
change the whole scenario of the laser plasma in- 
teraction [20 and, in particular, they can mod- 
ify the charged particle acceleration by elec- 
trostatic waves [211 El] and the injection [23] . 
The magnetic field with required for the injec- 
tion enhancement symmetry can be produced 
inside capillary plasma targets used for laser 
pulse guiding [24 . Imposed external homoge- 
neous magnetic fields can significantly improve 
the quality and stability of the laser wakefield 
accelerated electrons [25 . In the case of electro- 
static waves their frequency lies in the vicinity 
of the upper hybrid resonance [2 . Upper hybrid 
wake wave generation by laser light has been in- 
vestigated in Refs. [26]. Their properties have 
also attracted attention with regard to the ac- 
celeration of ultra high energy cosmic rays in 
the astrophysical environment [27j. 

When an extraordinary electromagnetic wave 
(whose electric field is perpendicular to the ex- 
ternal magnetic field, E ± Bq) of small ampli- 
tude propagates in the direction transverse to 
the magnetic field, k ± Bq, its frequency, cj, 
and wave number, /c, are related by the disper- 
sion equation 



Here ujpe = \/^i^ne^ /rUe and ujbc = eBo/rrieC 
are the Langmuir frequency and the electron 
Larmor frequency, respectively, with n, ^o, e, 
rrie and c being the plasma density, the ex- 
ternal magnetic field strength, electron charge 
and mass, and speed of light in vacuum. The 
ions are assumed to be at rest. When the mag- 
netic field vanishes, the dispersion equation (IT]) 
describes two independent waves: the electro- 
static longitudinal (k||E) Langmuir wave with 
the frequency equal to uj = uOp^i and the elec- 
tromagnetic transverse (k ± E) wave, for which 
the frequency and wave number are related as 

UJ = Jk^c^ + cj^g. In a plasma with a finite ex- 
ternal magnetic field, 5o 7^ 0, the extraordinary 
electromagnetic wave has both components of 
the electric field parallel and perpendicular to 
the magnetic field. 

In order to illustrate the extraordinary elec- 
tromagnetic wave dispersion properties, in Fig. 
[1] we plot the square of the refraction index, 
7V^ = k^(? jijp' , equal to 



N^ 



-uj^{2uj1^^ujIJ^uj:^ 



UJ- 



(CJ2 



■^^ 



■ut 



(2) 



as a function of the wave frequency for relatively 
large and small magnetic field. As we see the 
refraction index tends to infinity at the upper 
hybrid resonance, where the wave frequency is 
equal to 



^UH 



'Be' 



(3) 




FIG. 1. Square of the refraction index, N^ — 
k^c? /(jo'^ ^ as a function of the wave frequency nor- 



malized by ujp 
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In the vicinity of the upper hybrid resonance, 
which corresponds to the limit k ^ ojuh/c^ the 



wave is electrostatic, i.e. one can neglect the 
displacement current effects. 



Py 



-VxBq, 
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II. PLANAR UPPER HYBRID WAVE 

A. Relativistic upper hybrid oscillations 

In order to compare the properties of a non- 
hnear upper hybrid wave in the planar and axial 
symmetric geometries we at first consider a pla- 
nar electrostatic wave propagating perpendicu- 
lar to a constant homogeneous magnetic field 
directed along the z-axis, B = BqBz^ when the 
functions under consideration depend on time 
and the coordinate x. The electron hydrody- 
namics equations and Maxwell's equations yield 
for the electron density n, electron momentum 
p = Px^x -^Py^y a^d the electric field E^: 



dtn ^ dx{nvx) = 0, 



(4) 



(11) 



The solution to the continuity equation for 
the electron density reads: n(xo, t) = no/\dxQx\. 
Integration of Eqs. (10) and (11) yields Py = 
{eBo/c)^ + Fy{xo) and Ex = ^nrioe^ + ^x{xo), 
which shows that Py and Ex are related to each 
other via 



CPy Ex 

eBo Airnoe 



cPy jxo) ^ E,(xo) 
eBo Airnoe 



(12) 



Here the functions ^y{xo) and Ex{xo) are de- 
termined by the initial conditions. Assuming 
^y{xo) = and Ea^(xo) = and substituting 
them in the r.h.s. of Eq. ([9| we find that 
this equation can be written in the Hamiltonian 
form with the Hamilton function 

m.Px) = 



dtPx + VxdxPx = -eEx VyBo, (5) 



dtPy + VxdxPy = -VxBo, (6) 



dtEx + VxdxEx = ^TTTioeVx, 



(7) 
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-p'^ and V. 
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The nonrelativistic 



where Vx = cpx/Jm 

cpy I ^mlc^^pl^pl. 

limit has been analyzed in Ref. [28 . 

Changing from Euler coordinates x, t to La- 
grange variables xq , t we obtain the relationship 
between x and xq: 



X = Xq +^(xo,t) 



(8) 



with ^(xo, t) being the displacement of the elec- 
tron fiuid element from its initial position, xq. 
The electron velocity is equal io Vx = i- Here 
and below a dot denotes a differentiation with 
the respect to time in the Lagrange variables. 

In the Lagrange coordinates the system of 
equations ([s] - [t]) takes the form 



Px = -eEx 



-VyBo 



(9) 



y^ 



{eBon' 



pIc^ 



27rnoe^^^ 



(13) 



Since the Hamilton function (13) does not 
depend explicitly on time, the conservation 
oi l-L{^^Px) = h{xo) gives a relationship be- 
tween the electron momentum, px^ and the dis- 
placement, ^ along the trajectory determined 
by h{xo). The maximal electron momentum, 
Px,m = mec/T^^^-^, with jrn = h/nieC^, and 
the displacement amplitude, ^rm are related to 
each other as 



V2 



+ 7r) 



• 27m£| 



1. 



Here and below 



SB 



OJb, 



(14) 



(15) 



pe 



The value of /i(xo) determines the trajectory 
of Hamiltonian system ([9|. Its phase portrait, 
which coincide here with the contours of con- 
stant values of the Hamiltonian (13), is pre- 
sented in Fig. |2j 




FIG. 2. Contours of constant values of the Hamil- 



tonian (13) for £b = 0.75, uope — 1 and h — 
9,18,27,36,45. 



B. Cold Avave breaking limit for an upper 

hybrid wave travelling with constant 

velocity 

We assume that the wave propagates with a 
constant phase velocity Vph — c^ph (the group 
velocity of the upper hybrid wave in a cold 
plasma vanishes) in a cold plasma with immo- 
bile ions and homogeneous density no- In this 
case all dependent variables are functions of the 
variable X = 
obtain 



cPpht. From Eqs. i4\ m we 



no 



PphJrnl 



-vl 



-pI 



Px 



Py 



Ex Jml& ^pI^pI^ BoPy 



Pph^mlc^ ^pI^pI - p^ 
' BoPx 



(16) 



, (17) 



' f^ph ^Jmlc^ ^pI^pI-Px 



E'^ = Airnoe 



Pphyjrnlc^ ^pI^pI - Px 



(18) 



. (19) 



Here and below a prime denotes a differentiation 
with respect to the variable X. From Eqs. (18) 



and (19) we see that the relationship between 



the y-component of the electron momentum and 
the electric field: 



i^x _ PyC 

47rnoe cBq 



= constant. 



(20) 



If the constant on the r.h.s. of Eq. (20) van 



ishes, equations (17) and (19) can be cast in the 
form 



Pph Jmlc^ -^pI^pI-Px 
— / ^= ^ 



Px 



Pp 



-pI 



-pI 



(21) 



Px 



When the x-component of the electron ve- 
locity becomes equal to the phase velocity of 
the wake wave, i.e. Vx = cj3ph^ the denomina- 



tor on the r.h.s. of Eq. (21) becomes equal to 



zero. This singularity corresponds to the wave 
breaking (for details see Refs. [12] [29]) when 
the electron density and the gradients of all the 
functions tend to infinity. The singularity is 
reached at the maximum of the e lectron mo- 
mentum, pi)r — rUecPph/ \ h^ — P'iy^. Here it is 

taken into account that due to symmetry the 
^/-component of the electron momentum van- 
ishes at the breaking point. Considering the 
wave structure in the vicinity of the singular- 
ity at X = Xi)r we can easily find that in 
the nonrelativistic limit, when jSph <C 1 and 
pi)r ^ T^e^ph^ the electron momentum depends 
on 5X = X — Xiyr as 



^^e^vp ^^e^vp 



Px = 

\ Vph J ' 

In the ultrarelativistic case, when jph - 
» 1, the x— component of the elec 



Vv^ 



■/?.\ 



tron momentum is 



/g \l/3 

Pbr - rrieCj'^r ( ^f^P^ ) {kuH^Xf . (23) 



Here 



k 



uo 



UH 



^pe 



Be 

Iph ' 



(24) 



The value of 27r /kuH is an order of magnitude 
within the relativistic nonhnear upper hybrid 
wave length. 

In Fig. [3] we show the profile of the break- 
ing upper hybrid wave obtained from numerical 



integration of Eqs. (17 - 19) for the phase ve- 
locity of the wave equal to cPph = 0.95c and a 
magnetic field corresponding to ujBe/^pe = — 2. 
It is clearly seen that at the wave breaking the 
wave profile becomes singular with the gradi- 
ents of all the functions approaching infinity. 
The y-component of the electron momentum 
changes direction at X = Xi^r- The electron 



density according to Eq. (16) tends to infinity 
as n (X (5X~^/^. 




FIG. 3. Profile of the breaking upper hybrid wave 
for uJBe/o^pe — —2 and Pph — 0.95. 



III. CYLINDRICAL UPPER HYBRID 
WAVE 

A. The Lagrangian form of the equations 

describing cylindrical upper hybrid 

oscillations 

Here we consider axisymmetric cylindrical 
oscillations of the electron component of the 



plasma assuming the ions with homogeneous 
density no to be motionless, the homogeneous 
constant magnetic field is directed along the z- 
axis, B = Bq^z^ and all the functions depend 
only on time t and the spatial coordinate r. In 
this case from the equations of the electron hy- 
drodynamics and Maxwell's equations we ob- 
tain the following system: 



dtn -\ — drirnvr) = 0, 
r 



(25) 



dtPr + VrdrPr 



PifVjf 



-eEr--v^Bo, (26) 



dtP^ + VrdrP^ 



-^ = -VrBo, 



(27) 



dtEr 



- Vr-drirEr) = ^TiriQeVr- (28) 

r 



Here n is the electron density, Pr and p^ are 
the radial and azimuthal components of the 
electron momentum, and Er is the radial com- 
ponent of the electric field. The radial and 
azimuthal components of the electron velocity 

cpr/^Jrril 



are equal to Vr 



Pr 



p'^ and 



Vy, = cp^/ Jml 



-Pr^p'^^ respectively. 
One of the most distinct difference between 
planar and cylindrical geometry is the appear- 
ance of the centrifugal force in the l.h.s. of Eq. 
(26). The azimuthal momentum changes due to 



interaction of the radial motion with the mag- 
netic field and vice versa. The effect of the cen- 
trifugal force becomes stronger the closer the 
radial flow of the electron component converges 
thus preventing the flow collapse on the axis. 

In order to solve Eqs. ([26] - [28| ) it is conve- 
nient to change from Euler coordinates, r, t, to 
Lagrange coordinates, ro,t, [17]. The Euler and 
Lagrange variables are related to each other as 



r = ro + p(ro,t). 



(29) 



where p(ro,t) is a displacement of the elec- 
tron fluid elements from its initial position, i.e. 

p(ro,0) = 0. 



In the Lagrange coordinates the continuity 
equation (l25| has a solution 



n(ro,t) = 



no 



\J{ro.tW 



(30) 



where J(ro,t) = (ro + p)(l + droP)/ro is the 
Jacobian of the transformation from the Euler 
to Lagrange variables. Equations ([26]- 28) can 
be written in the form 



Pr 



PifVif 



-eEr Vn,Bo, 



(31) 



H r (f) — TTTioe r -\- e 

2c 



In—. (37) 

^0 



Since the Lagrange function ( 37 ) does not de 



pend on time explicitly, the conservation of the 
Jacobi integral, 



dC 



Qj ->C = ^(ro), 



(38) 



Qj 



where Qj = (r, (p) and Qj = (r, (^), results in the 
energy conservation 



P^ 



Er 



PifVr 



ErVr 



-VrBo, 



AnrioeVr. 



(32) 



(33) 



Taking into account that the radial and az- 
imuthal components of the electron velocity are 
equal to Vr = p and v^ = np with cp being the 
azimuthal angle, we find solutions to Eq. (32) 
for the azimuthal momentum 



^P^ = ^^or 



:'(^0) 



and to Eq. (33) for the electric field 
rEr = 27rnoer^ - 



(34) 



(35) 



These expressions are the consequences of the 
conservation of the azimuthal component of the 
canonical momentum and of the charge conser- 
vation, respectively. Here functions M^{ro) and 
Q(ro) are determined by the initial conditions. 
Similar to the planar geometry case, when the 
electric field and the ^/-component of the elec- 
tron momentum are related to each other via 



Eq. (20), in cylindrical geometry we have 



Erp 

47rnoe 



P^c \ 



Q(rQ) 



.{ro) 



eBo 



(36) 



Equations ( |31|) a nd (|32| with the electric field 
given by Eq. (35) can be written in the La- 



grangian form with the Lagrange function 



£ 



-rripC 



+ r 



2^2 



c 2 2 

t = —irnoe r 



ln-+^(ro). (39) 

^0 



On the l.h.s. of this equation the energy is equal 
to 



£ 



pI 



-pi- 



rn 



In the r.h.s. the function ^(ro) is determined 
by the initial conditions. Using expressions (34) 



and ( 39 ) it is easy to find the radial component 



of the electron momentum pr as a function of 
the coordinate r. Then with the equation for 
Vr, which is 



2-Pr 



8' 



(41) 



we solve by quadratures Eqs. (31 -33), because 
the r.h.s. of Eq. ( [41] ) explicitly depends on the 
coordinate r only. Its dependence on the La- 
grange coordinate is parametric. For a known 
time dependence of the radial coordinate, r(t), 
integration on time of the equation for the az- 
imuthal angle 



2P^ 



(42) 



yields the trajectory of the element of the elec- 
tron fiuid. Its projection to the r, (p plane is 
given by 



^■ 



Pipdr 



constant. 



rpr 



(43) 



In the case when the initial value of the 
electric field and azimuthal momentum vanish, 
the functions Q(ro) and M(^(ro) are equal to 



Q(ro) = -2^noer2 andM^(ro) = -{eBQ/2c)rl, 
respectively. If the radial monienturQ initial 
value equals Pr,mi from Eq. (39) it follows 



^0 



mld^ + 



eBc, 



-^vic 



2^2 



mld^ ^Plrn^^ 



(44) 
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Fig. [4] shows the energy iso-contours plotted 
for expression ( [44| ) with ujbc = 2 and ujpe = 1 
As we see the radial coordinate cannot vanish 




FIG. 4. Iso-energy contours for expression (44) for 

ujBe — 2 and ujpe — 1 and /i = 9, 18, 27, 36, 45. 

due the logarithmic divergence of the electro- 
static potential and the ex 1/r divergence of the 
"centrifugal energy" at r ^ 0. The minimal 
value of the radius can be estimated to be ap- 
proximately equal to 



e^o 



(45) 



Fig. [5] presents the typical electron trajec- 
tory in the {x = r cos cp^y = rsimp) plane for 



OJBe/^: 



1.5, pr{0) = -0.5, p^{0) 



-0.25, and 



pe 



0.25. It is clearly seen that the 
trajectory is localized in the region r > Tmin, 
which, in the case under consideration, is about 
0.31. 





2 


y 






P^ 


\j\)( / 


%, 


-2(^^^-l 




^^2 








^^" 




/ /(Y) 


w 



FIG. 5. Electron trajectory in the [x — r cos Lp^y 
rsimp) plane for ro = 1.5, Pr(0) = —0.5, Pc^(O) 
—0.25, and uJBe/ojpe = 0.25. 



B. On the cavity closing in a Avake left 
behind the laser pulse 

As is well known, the cavity, the positively 
charged region void of electrons, is formed be- 
cause almost all the electrons in the way of a 
sufficiently intense laser pulse are pushed aside 
[m [TQ . The unperturbed electrons on the side 
are attracted by the positive charge. They close 
the cavity at the finite distance behind the laser 
pulse. A non-zero axial magnetic field prevents 
the cavity from closing as shown in the previ- 
ous section. The process of cavity closing is ex- 
plained in Fig. [6Jin the case of vanishing mag- 
netic field (Fig. 16] a) and for finite magnetic field 
(Fig. [6]b) when it results in the appearance of 
a hole at the rear of the cavity. 

In order to estimate the radius of the hole at 
the rear wall of the cavity we choose the pa- 
rameters of the converging cylindrical wave as 
follows. The functions M(^(ro) and Q(ro) deter- 



mining the azimuthal momentum in Eq. (34) 



and the radial component of the electric field in 




FIG. 6. The cavity formed behind the laser pulse in 
a collisionless plasma. The walls formed by unper- 
turbed electrons converge at the rear of the cavity 
(a). In the presence of a longitudinal magnetic field, 
the electrons miss the axis, which results in the ap- 
pearance of a hole at the rear of the cavity (b) . 



(|35]), respectively, are assumed to be equal 
^(^(^o) = and 



[ro) 



fo, 

\27 



27rnoe(r2-rg), 



ro < r^; 
ro > Tc. 



(46) 



Here Vc is the initial cavity radius. Using this 
expression we obtain that at t = the electric 
field is a linear function of radius for ro < Vc 
and it is inversely proportional to ro at ro > Vc 
(see Fig. [7|), 



Er{ro) = — 27rnoe 



^0, 



ro < r^; 



r^/ro, ro > r^, 



(47) 



in accordance with the radial electric field pro- 
file found with the 3D PIC simulations (see Fig. 
2 d,e in Ref. ^16j). Considering motion of the 




FIG. 7. Electric field as a function of the Lagrange 
coordinate, ro. 

electron fluid element with ro = Vc and the ini- 
tial value of the radial momentum equal to zero 



we obtain from Eq. (39) 




'_c_ 

r 



4c2 



{rl-r') 



(48) 



According to Eq. ( |48[ ) the radial component of 
the electron momentum vanishes at r = re and 
at r = rmin, which is given by equation 




4c2 



pe / 2 



(49) 



In Fig. [8] we plot dependence on the mag- 
netic field {(jJBe/^pe) of the minimum radius, 
^min, normalized by 2de = c/ujpe^ for different 
values of the cavity radius, Tc. In the limit of 




«../«,. 



FIG. 8. Minimum radius, rmin, normalized by de = 
c/uope v.s. the parameter ( 15 ) characterizing the 
magnetic field amplitude, for the cavity radius, re = 
1,2,3,4,5,6,7 

weak magnetic field, when ooBe/^pe "C 1, the 
minimal radius, rmin is linearly proportional to 
ojBe^ being approximately equal to 



4 CUJBe 

V2 Ujle 



(50) 



i.e. it does not depend on the initial cavity 
radius, Tc- For B ^ the hole radius tends 
to zero. In the limit of strong magnetic field, 
^Be/^pe ^ 1, we have Tmin ^ ^c, i-G. it wcakly 
depends on the magnetic field amplitude. 

Fig. |9] presents the typical electron trajec- 
tory in the {x = r cos (p^y = rsimp) plane for 
Tm = ro = 1, Pr{0) = 0, p^{0) = 0, and 
^Be/^pe = 0.5. It is clearly seen that the tra- 
jectory is localized in the region r > rmin- 




FIG. 9. Electron trajectory in the (x = r cos (p,y = 
rsinc^) plane for rm = To = 1, Pr{0) = 0, P(^(0) = 0, 
and (jJBe/(^pe = 0.5. 

By using values of the Langmui r and Larm or 
frequencies, ujpe = 5.642 x 10^ ^/ no {cm~^)s~^ 



and UBe = 1.759 x 10^1B{T)s-^ 
the minimal radius for no = 10^^cm~ 
lOT to be Tmin = 4.677 x 10-^cm. 



we estimate 



and B 



C. Betratron oscillations of relativist ic 

electrons inside the Avake Avave cavity in the 

presence of the axial magnetic field 

The effects of axial magnetic field on the 
LWFA has been addressed in Refs. [30] and [3T] 
where the plasma densification through mag- 
netic compression has been studied in order to 
overcome dephasing between the wake field and 
accelerated electrons and the magnetic field ef- 
fects on the electron injection have been anal- 
ysed, respectively. Here we pay attention to the 
electron trajectory rotation in the azimuthal di- 



rection resulting in the azimuthal patterns of 
the LWFA accelerated electron bunches which 
has been observed in the experiments reported 
in Ref. [25l. The betratron oscillations in the 
laser wake field themselves provide a promis- 
ing source of ultra short x-ray beam generation 
[32 . The magnetic field can be used for addi- 
tional manipulation of the betatron oscillations 
and for controlling the properties of the radia- 
tion emitted. 

Inside the cavity formed in a plasma in a wake 
behind an ultra short laser pulse the electric 
field has a radial component Er = 2iTner. Fol- 
lowing betatron oscillation theory [33 we as- 
sume that the longitudinal (axial) component 
of the electron momentum, p^i is given. Choos- 
ing as a function of time 



P.(i)=Pn^(iacc -*')/*' 



(51) 



we describe the electron acceleration in the 
wake field [34 . During the acceleration time 
^occ ^ {^/^pe)lph ^^^ electron longitudinal mo- 
mentum changes from zero to Pm ~ rrieC^p^ 
[5l. Here the gamma factor associated with 
the phase velocity of the wake wave is equal 
to jph = ^/uJo/uJpe. The equations of trans- 
verse electron motion are identical to Eqs. (31) 



and (32) where the radial electric field is given 
by Er = 27Tner. The electron velocity is v = 
{Px/l^Pr/l^Pif/j) with the relativistic gamma 



pI' 



factor 7 = w 1 + p^ ^ ^2 

In the case when the magnetic field vanishes, 
the electron trajectory is described by a linear 
combination of associated Legendre functions 
[34]. If the electron is injected into the wake 
wave with zero azimuthal momentum its trajec- 
tory lies in the same plane. The axial magnetic 
field causes the electron trajectory to rotate as 



seen in Fig. [10) In this figure we plot the elec- 
tron trajectory obtained via numerical integra- 
tion of the electron equations of motion for a 
normalized magnetic field B = 0.25, initial co- 
ordinates, ro = l^ (po = 7r/2, xq = —800, initial 
momentum p^^ = 0, p^ = 0, P(^=o, acceleration 
time tacc = 800 and maximal longitudinal mo- 
mentum equal to to Pm = 50, respectively. 
We distinctly see the betatron oscillations 
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FIG. 10. Betatron oscillations for B = 0.25, ro = 1, 
(fo = 7r/2, xo = -800, px = 0, pr = 0, p^^o, 
tacc — 800, and pm — 50. The electron trajec- 
tory (black curve) is confined in a slowly rotating 
plane. The trajectory projections on the bottom 
of the X, y plane (green curve) and on the right x, z 
plane (blue curve) show slowly varying amplitudes 
and periods of the betaron oscillations. In the up- 
per X, y plane we plot the longitudinal momentum 
Px as a function of the x-coordinate (yellow). The 
trajectory projection on the r.h.s. ?/, z plane (red) 
shows a "star"-like pattern. 




FIG. 11. Images of the electron beam generated 
with B = 0.2T, gas density riHe = 4 x lO^^cm^, and 
a laser power of 12 TW. The laser polarization is 
in the vertical direction, a) Four-ray star, b) Spiral 
pattern. 



of the electron whose amplitude and period 
changes according to dependence of the beta- 
tron frequency on the electron energy, cj^ = 
cjpe/T, and by virtue of the adiabatic invari- 
ant conservation fperp/^6 = constant [34j. Here 
Eperp is the kinetic energy of the transverse mo- 



tion. The electron trajectory (black curve) is 
confined on a slowly rotating plane. The tra- 
jectory projections on the bottom of the x, y 
plane (green curve) and in the right x, z plane 
(blue curve) show slowly varying amplitudes 
and the period of the betaron oscillations. On 
the upper x, y plane we plot the longitudinal 
momentum px as a function of the x-coordinate 
(yellow). The trajectory projection on the r.h.s. 
y^ z plane (red) shows a "star" -like pattern. We 
note an analogy of this pattern with the images 
of the electron beam generated in the exper- 
iments [25] on LWFA with an axial magnetic 
field. The "star-like" pattern should also be re- 
vealed in the betatron emission. In Fig. pTjwe 
present images of the LWFA electron beam ob- 
served in [25 for the axial magnetic field equal 
to 0.2 T. We see a "Four-Ray Star" image in 
Fig. fll] a) and a spiral pattern in Fig. 11 b), 
which clearly indicates the magnetic field effect. 
When the magnetic field vanishes a typical 
electron beam image has an elliptic form, as in 
Fig. 12 The electron trajectory (black curve) 



has the form of a spiral with a slowly chang- 
ing amplitude and period. The trajectory pro- 
jections on the bottom of the x, y plane(green 
curve) and on the right x, z plane (blue curve) 
show slowly varying amplitudes and the period 
of the betaron oscillations. On the upper x, y 
plane we plot the longitudinal momentum p^ 
as a function of the x-coordinate (yellow). The 
trajectory projection on the r.h.s. y, z plane 
(red) shows an "ellipse" -like pattern. Such el- 
liptic structures have been discussed in detail 
and have been detected in experiments [35]. 



IV. COMPUTER SIMULATION 
RESULTS 

We note here that the above analyzed prop- 
erties of nonlinear plasma waves in magnetized 
plasmas cannot be revealed in a 2D planar (x, y) 
symmetry configuration due to the vanishing of 
the centrifugal force. In order to reveal the 3D 
effects in the configuration typical for the ex- 
periments in high power ultrashort pulse laser 
interaction with underdense plasmas we per- 



11 




FIG. 12. Betatron oscillations for B = 0, ro = 1, 
ifQ = 7r/2, xq = -800, px = 0, pr = 0, Pc^=i, 
^acc = 800, and pm — 50. The electron trajectory 
(black curve) has the form of a spiral with slowly 
changing amplitude and period. The trajectory 
projections on the bottom of the x, y plane(green 
curve) and on the right x, z plane (blue curve) show 
slowly varying amplitudes and the period of the be- 
taron oscillations. On the upper x^y plane we plot 
the longitudinal momentum px as a function of the 
x-coordinate (yellow). The trajectory projection in 
the r.h.s. y^z plane (red) shows an "ellipse" -like 
pattern. 



(a,d), a close-up of the isosurface correspond- 
ing to ne/ucr = 0.008 is shown in the location 
of the wave-breaking (b,e); the cross section is 
about X = 60. 9A averaged over a half of the laser 
wavelength (c,f). In Fig. [l3] (a,b,c) we see the 
cavity formed in the electron density with thin 
high density walls. The rear side of the cavity is 
closed resulting in the electron bunch injection 
along the axis. The electron density has a local 
maximum on the axis (Fig. fl3](b,c)). The pres- 
ence of the static magnetic field prevents the 
cavity from closing (Fig. IT3](a,b)) in accordance 
with the above formulated theory of the cylin- 
drical Upper Hybrid plasma wave. The elec- 
tron bunch injected along the cavity axis from 
its rear side is ring-shaped with no electrons on 
the axis. The electron density cross section re- 
veals the "Four-Ray Star" pattern which has 
been observed in the experiments [2F. 

In the region relatively far from the axis 
where the magnetic field is relatively small, Fig. 
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(a,d), the electron density distribution is not 
affected by the magnetic field showing identical 
shapes of the cavity walls and bow wave. 



CONCLUSION 



formed three dimensional Particle in Cell simu- 
lations with the code REMP [36 whose results 
are presented in Fig. [13] The Gaussian lin- 
early polarized (along z-axis) laser pulse with 
the amplitude of ao = 6 and FWHM dimensions 
of 8A X 20A X 20A propagates in homogeneous 
plasma with the density of ng = 0.0009ncr • Ions 
are assumed to be immobile. The simulation 
box size is 80A x 72A x 72A, the mesh sizes are 
dx = A/16, dy = dz = A/8, and the total num- 
ber of quasi-particles is 6.3 x 10^. 

Fig. [13] (a,b,c) presents the electron density 
in the wake of the laser pulse propagating from 
left to right in a homogeneous plasma, when 
no static magnetic field is imposed. In Fig. 
[13] (d,e,f ) we plot the electron density distribu- 
tion in the presence of a homogeneous perma- 
nent longitudinal magnetic field of 10 T. Half 
of the box is removed to reveal the interior 



The effect of the axial magnetic field on the 
cylindrical electrostatic wave dynamics leads to 
the precession of the electron trajectory in the 
(x = r cos (p^y = r sin cp) plane. Conservation of 
the generalized momentum leads to a non-zero 
angular momentum of the electron component. 
As a result the electrons do not reach the axis 
remaining confined in the region r > rmin thus 
preventing the rear wall of the cavity in the 
wake behind ultra-short laser pulse from clos- 
ing. In the limit of weak magnetic field, when 
^Be/^pe <^ I7 the minimal radius, r^i^ is lin- 



early proportional to uJBe^ as given by Eq. (50). 
The magnetic field causes a rotation of the fast 
electron trajectory, which is seen in a typical 
"Four- Ray Star" pattern of the electron density 
image observed in the experiments [25 (see also 
Fig. 11 and in the 3D PIC simulations in Fig. 



13). 
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FIG. 13. Electron density in the wake of the laser pulse propagating from left to right in a homogeneous 
plasma with no static magnetic field initially imposed (a,b,c) and in the presence of a longitudinal magnetic 
field of 10 T (d,e,f). Half of the box is removed to reveal the interior (a,d), a close-up of the isosurface 
corresponding to ne/ricr = 0.008 is shown in the location of the wave-breaking (b,e); the cross section is 
about X = 60. 9A averaged over a half of the laser wavelength (c,f). 
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